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We theoretically study Hawking radiation processes from an analog acoustic black hole in a 
flowing superfluid of exciton-polaritons in a one-dimensional semiconductor microcavity. Polaritons 
are coherently injected into the microcavity by a laser pump with a suitably tailored spot profile. 
An event horizon with a large analog surface gravity is created by inserting a defect in the polariton 
flow along the cavity plane. Experimentally observable signatures of the analog Hawking radiation 
are identified in the scattering of phonon wavepackets off the horizon, as well as in the spatial 
correlation pattern of quantum fluctuations of the polariton density. The potential of these table-top 
optical systems as analog models of gravitational physics is quantitatively confirmed by numerical 
calculations using realistic parameters for state-of-the-art devices. 
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I. INTRODUCTION 

The quantum mechanical properties of the vacuum 
state in quantum field theories are presently attracting 
the attention of researchers from very different communi- 
ties, from astrophysics and gravitation, to quantum op- 
tics and condensed matter physics. A central idea of all 
these studies is the possibility of converting zero-point 
quantum fluctuations into observable quantum vacuum 
radiation by some spatial and/or temporal dependence 
of the background over which the quantum field is prop- 
agating^. A number of fundamental physical effects be- 
long to this category, from the dynamical Casimir effect 
when the boundary conditions of the quantum field are 
rapidly varied in time^ ^, to the Hawking radiation when 
it prop aga tes on a curved space-time showing a black hole 
horizon^. 

In the wake of the pioneering work by Unruh^, re- 
searchers have started investigating condensed-matter 
systems where the propagation of some low-energy ex- 
citation field follows an effectively curved space-time ge- 
ometry. In suitable configurations showing a black hole 
horizon for this low-energy excitation, a quantum vac- 
uum emission is expected to appear via a mechanism 
analogous to Hawking radiation. Among the many sys- 
tems that have been investigated in this perspective^, 
low-temperature superfiuids and nonlinear optical sys- 
tems are nowadays considered the most promising ones. 

Dilute superfiuids such as Bose-Einstein condensates 
of ultracold atoms join a very simple spectrum of ele- 
mentary excitations with sonic dispersion, with the pos- 
sibility of pushing the sensitivity of measurements down 
to the quantum limit. Their potential as analog mod- 
els was first proposed in Refs. [TQHT2| and later confirmed 
by ah initio numerical simulations of the condensate dy- 
namics'^. Even if acoustic black hole configurations have 
been experimentally realizecJ^, no evidence of Hawking 
radiation has been reported yet. 



Analog models based on nonlinear optical systems were 
pioneered in Ref. [151 a strong pulse of light propagating 
in a Kerr nonlinear medium can be used to generate a 
moving spatial interface separating regions of different 
light velocity. For suitably chosen parameters, the inter- 
face then behaves as a horizon that can be crossed by 
light in one direction only, which should result in Hawk- 
ing radiation being emitted. Some experimental claim 
in this direction was recently reported'^ and has raised 
a number of interesting questions concerning the inter- 
pretation of the observed radiation in terms of Hawking 
processeJI^KI^. 

A completely new perspective to analog models based 
on nonlinear optical systems was opened by MarincP^, 
who first proposed the use of quantum fiuids of light 
to generate acoustic black hole configurations and, 
then, to look for the Hawking radiation of Bogoliubov 
phonons on top of the photon fiuid. In the following 
years, this idea has been pushed forward by a number 
of authors, who h ave considered different geometrical 
configurationP^^ and material systems, and very re- 
cently exciton-polariton fiuids in semiconductor micro- 
cavitie^. 

In the present work, we report a comprehensive study 
of analog Hawking radiation effects in quantum fiuid of 
light, specifically exploiting the quantum fiuid proper- 
ties of an exciton-polaritons condensate. Our theoretical 
model fully includes their intrinsically non-equilibrium 
nature^' and describes quantum fiuctuations of the po- 
lariton field within the so-called truncated- Wigner for- 
malism of degenerate quantum gaseP^'^. As a result, 
our numerical calculations are able to provide quanti- 
tative predictions for the observable quantities, and to 
point out clear and accessible signatures of analog Hawk- 
ing radiation in the emitted light from the cavity. Even 
if our discussion is mostly focused on the specific case 
of semiconductor microcavity device^^^ where po- 
lariton superfiuidity has been first demonstrated^, all 
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our conclusions straightforwardly extend to generic pla- 
nar cavity devices filled with or made of a Kerr nonlinear 
medium. 

The article is organized as follows. In Sec. |ll| we in- 
troduce the physical system under consideration and we 
summarize the theoretical tools that are used to describe 
it. In Sec. Ill we describe the original laser beam con- 



laser pump 



figuration that we propose to generate an analog black 
hole horizon with a large surface gravity. The observable 
consequences of Hawking mode-conversion processes are 



illustrated in the following sections. In Sec. [IV| we dis- 
cuss the scattering of a coherent phonon wavepacket on 
the horizon: the signature of classical Hawking processes 
is visible as an additional wavepacket emerging from the 
horizon on a negative-norm branch. In Sec.[V| we present 
the numerical evidence of the Hawking radiation due to 
the conversion of zero-point quantum fiuctuations into 
observable radiation, and we compare it with the the- 
oretical expectations. Conclusions are finally drawn in 
Sec. ED 



sub-sonic super-sonic 
flow horizon flow 



analog Hawking emission 



> 

CD 

E 



CD 
C 
LU 



1590 



1585 



1580 



1575 



II. PHYSICAL SYSTEM AND THEORETICAL 
MODEL 

We consider a polariton wire device where cavity pho- 
tons propagating along one dimension are resonantly cou- 
pled to the fundamental exciton transition in one or more 
quantum wells embedded in the cavity layer (in-plane po- 
larized heavy- hole to conduction band transition). Start- 
ing from a planar microcavity where light is confined 
along the growth axis z by a pair of distributed Bragg 
refiectors, a suitable etching procedure is used to later- 
ally pattern the device in the y direction, and create the 
ridge structure that is schematically shown in Fig. [TJa). 
Provided the losses are weak enough, the bosonic excita- 
tions that result from the strong coupling of the cavity 
photon with the quantum well exciton have the typical 
mixed light-matter nature of one-dimensional exciton- 
polaritons, free to propagate along the wire axis accord- 
ing to the dispersion law that is plotted in Fig. [ijb). 
A high quality experimental realization of this polariton 
wire concept was reported in Ref. 33 , from which we shall 
take the state-of-the art system parameters to be used in 
the numerical simulations. 



A. The quantum field Hamiltonian 

A theoretical description of the dynamics of this sys- 
tem can be developed in terms of the standard Hamil- 
tonian for the two coupled bosonic fields describing the 
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Figure 1: (Color online) Upper panel: sketch of the polari- 
ton wire device under consideration. Lower panel: energy- 
dispersion of the upper and lower polariton branches as a func- 
tion of the wavevector kx along the wire axis. The thin dashed 
and dash-dotted lines indicate the bare cavity photon and ex- 
citon dispersions, respectively. The thick dashed line indicates 
the parabolic approximation of the lower polariton dispersion. 
System parameters are inspired from Ref.l33l huj^^ — hwi^^ — 
1580 meV, hTtR - 7.5 meV, nic - 1.2 x 10~^mo (with mo the 
free electron mass), ^ rric. 



quantum well exciton and the cavity photorP^ 

+ dx \IjI{x)^1{x)iIj^{x)^^{x) 

+ /dx«E(.,,W,.M + h.c.. (1) 

where x is the longitudinal coordinate along the wire axis 
and the two-component operator vector ^ = (V^x,'^c)^ 
summarizes the quantum fields describing the exciton 
?/^x(^) and photon V^c(^) fields. Each of them sat- 
isfies one-dimensional bosonic commutation rules, e.g. 
[Vii(x), Vi](x')] = 5i^j5{x - x'), with i,j = {x,c}. 
Throughout the whole paper, we restrict ourselves to the 
case where the system is pumped on a single spin state. 
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so that the spin degrees of freedom can be neglected in 
the theoretical model. 

The single-particle Hamiltonian Hq describing the evo- 
lution of the non-interacting exciton and cavity photon 
fields has the simple representation 

TT _ ( hu^{-idx) ^ V^{x) hnn \ . . 

in terms of the bare cavity photon and exciton dispersion 
law ujc,^{kx) ^ cjc^x + h^kx^ / {2mc^^) in a spatially homo- 
geneous one-dimensional system. The rest energy uj^^ 
and the effective mass mc of the cavity photon are deter- 
mined by the spatial confinement along the z and y axis 
by the DBR mirrors and by the total internal reflection 
at the etching interfaces, respectively. In practical cal- 
culations the exciton dispersion can be safely neglected, 
^^{kx) = as the exciton mass is orders of mag- 
nitude larger than the photon mass. The strength of 
the light-matter coupling is quantified by the Rabi fre- 
quency Qr, proportional to the square root of the oscil- 
lator strength per unit area of the excitonic transition^^. 

In spatially homogeneous systems, for which the exter- 
nal potentials acting on the cavity photon and the quan- 
tum well exciton vanish {Vx,c{x) = 0), diagonalization of 
the single particle Hamiltonian Hq as a function of wave 
vector kx leads to the well-known upper and lower polari- 
ton branches of dispersions 0Jup{kx) and ojLp{kx)- As it 
is shown in Fig. [ijb) the bottom part of the lower polari- 
ton dispersion is accurately approximated by a parabolic 
form 



^Lp{kx) 
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where the uj^^p and rriLP parameters have the clear phys- 
ical meaning of a rest energy and an effective mass for 
the (lower) polariton. In the resonant case uo^^ = uji^\ 
one has the simple relation tulp — 2mc. 

Several methods to generate external potentials 14,c(^) 
with almost arbitrary spatial shapes along the x direc- 
tion have been experimentally demonstrated in the last 
years, from the in-plane patterning of the microcavity 
layei!^, to the all-optical potential created by a strong 
pump laser with opposite polarizatiorP^, to the applica- 
tion of a mechanical stress to the device^. 

Exciton-exciton interaction, due to the Coulomb inter- 
actions between the electron and holes forming the exci- 
tons, are the main source of optical nonlinearity in this 
problem. Experimental measurements^^ report a value 
on the order of %2D ^ 2.5 /ieV • /im^, which translates 
into a reduced one-dimensional interaction parameter to 
be used in the Hamiltonian (|l|) ^ = g2D J^^y\^x{y)\^j 
where the integral is over trie transverse wire width 
(— I/w/2 < y < I/w/2), and the forth power of the ex- 
citon wave function comes from the Kerr-type nature 
of the nonlinearity. Assuming the transverse confine- 
ment of the exciton envelope function to be described 



by ^x{y) = \/2/Lw cos {ny/L^)^ a good estimate can 
be given to be ^ = 6^2d/^w Taking into account the 
exciton weight in the polariton, this translates in the res- 
onant uj^^ = cjc^^ case into a one-dimensional polariton- 
polariton interaction constant qlp g/A. 

Polaritons are coherently injected into the system by 
an external laser drive coupled to the cavity photon via 
the non-perfect reflectivity of the DBR mirrors. In the 
theoretical model, such processes are described by the 
terms in the last row of Eq. ([T]), where E{x^t) is the 
spatio-temporal profile of the coherent pump beam in- 
cluding the coupling coefficient proportional to the DBR 
transmission amplitude. Correspondingly, photons (ex- 
citons) are subjected to radiative losses at a rate 7c (7x) 
that have to be described at the level of the master equa- 
tion for the density matrix, p, for the coupled quantum 
fields 



dp 
dt 



ih 



[n,p] + c[p]. 



(4) 



Performing the standard Born-Markov approximatioiP^, 
the dissipation superoperator can be cast into the Lind- 
blad form 



2pi{x) p'lplix) - 2pl{x) 2pi{x) p 

- pipl{x)ipi{x) . (5) 



The resulting decay rate of polaritons is a weighted 
average of the exciton and photon ones 7x,c* in the 
vicinity of the bottom of the lower polariton branch of 
Fig.jljb), it has an almost momentum-independent value 
7LP ^ (7x + 7c)/2. 



B. The generalized Gross-Pitaevskii equation 

In a regime of weak exciton-exciton interactions, the 
dynamics of the system can be accurately captured by 
the mean-field approximation: the quantum fields ?/^x,c 
are replaced by the classical fields corresponding to their 
expectation values 



^x,c(^) = (V^x,c(^)) 



(6) 



and evolving according to the pair of nonlinear partial 
differential equations 



ih— 
dt 

Hq - 



(t)^{x,t) 

hg\(l)^{x,t)\^ 






hE{x,t) 



^{x,t) 



(7) 



which generalize to the non-equilibrium context of polari- 
tons the well-known Gross-Pitaevskii equation (GPE) of 
dilute Bose condensed gases^^. 
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wavevector kp are indicated here as Acp = ujc{kp) 
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Figure 2: (Color online) Upper panel: hysteresis loop in the 
polariton density vs. incident laser amplitude for a plane wave 
pump of wavevector kp = 0.2 /im~^, and frequency hwp = 
1574 meV tuned slightly above the lower polariton branch at 
(jOLp{kp) 1572.56 meV. Equal loss rate are assumed for the 
exciton and the photon fields, hj^^c = 0.02 meV. The value 
hg = 5 /ieV • /im of the nonlinearity corresponds to a Lw = 
3/im wide polariton wire (see text). Lower panel: spatial 
intensity profile of the coherent pump along the wire axis, 
with excitation parameters: HFi — 24.7487 meV • fim~ 
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hF2 1.1314meV•/im"^/^ kp 0.2/im 



1574 meV. 



and Axp = uj^{kp) - Up. 

In the quantum fluid language, this equation plays the 
role of the equation of state for the polariton fluid under a 
coherent pump in a plane wave form: while the polariton 
flow speed is controlled by the incident wavevector /cp, 
the density has a more complex dependence on the pump 
intensity and frequency ujp. A most remarkable 

example is illustrated in Fig. [2^: the choice of a pump 
frequency slightly above the lower polariton branch leads 
to an optical bistable behavior. 

For reasons that will soon become clear, we shall con- 
centrate our attention onto the case where the pump in- 
tensity is tuned in the close vicinity of the end-point of 
the upper branch of the hysteresis loop (indicated as F2 
in Fig. |2|i). In order to reach this working point in a 
given spatial region, a pump with a non-trivial spatial 
envelope, Fp{x)^ multiplying the plane wave can be used, 
as shown in the lower panel of Fig. [2|b). Upstream of 
the region of interest (i.e. for 10 fim < x < 50 fim the 
pump amplitude is tuned to a value Fi higher than the 
switch-on point of the hysteresis loop (i.e. the end point 
of the lower branch), so to switch the system to the up- 
per branch of the hysteresis loop^. In the long region of 
interest {50 jam < x < 350 /im), the pump amplitude is 
maintained at a spatially constant value F2, so to keep 
the system sits at the desired working point with the 
desired in-plane flow wavevector k^. 

The behavior of the density n{x) and the wavevector 
k{x) in the downstream region (350 /im < x) is deter- 
mined by a complex interplay of propagation and losses. 
First numerical and analytical studies of this regime have 
appeared in Refs. '42^and'43^ the spatial drop of the po- 
lariton density away from the pump results in an increase 
of the polariton speed under the effect of the interactions. 
The lower the losses, the slower this acceleration effect. 



In the simplest case of a spatially homogeneous system 
under a coherent pump in a plane wave form 



E{x^ t) = Fp ex.-p[i{kpX — ojpt)] 



(8) 



of frequency cjp, wavevector kp and amplitude Fp, the 
analog of the equation of state can be obtained by sim- 
ply injecting the same plane wave ansatz (j)x,c{x^t) = 
^J^^ exp[i{kpX — ujp)t] into the generalized GPE 
This leads to the simple form 



(9) 



for the steady-state excitonic density = as a 

function of the pump parameters, and to an analogous 
expression for the cavity-photon densit^p^^. The pump 
detuning from the cavity photon and exciton modes at 



C. Bogoliubov dispersion of collective excitations 

The dispersion of collective excitations on top of the 
steady state of the coherently pumped condensate is ob- 
tained as usual by linearizing the GPE around the steady 
state solution^^. A complete discussion of the many dif- 
ferent behaviors that are possible depending on the pump 
frequency and intensity can be found in Refs. I25|26l 
Here, we shall restrict our attention to the two cases of 
present interest, namely (i) a coherent pump whose in- 
tensity is tuned right at the turning point of the upper 
branch of the hysteresis loop, and (ii) the ballistic flow 
regime in the downstream region. 

If the polariton population is restricted to the lower 
polariton branch and we perform a parabolic approxima- 
tion of the branch bottom with effective mass m^p, the 
dispersion of the collective excitations on top of a con- 
densate of polariton density n^p in uniform flow at speed 
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vlp has, in both cases, the usual BogoUubov form 

l h{k,-kLPY f h{k,-kLPY ^ ^^^^ 

y 2mLP V 2mLP ) 

^VLp{k,-kLp)^'-^. (10) 

In the first (i) case, the local wavevector kLP is fixed to 
the in-plane component kp of the pump w avevect or. In 
the second (ii) case, hydrodynamic effectj^ ^^*^^^ make 
the local wavevector kLp{x) to slowly vary in space. In 
both cases, the local in-plane flow speed vlpJs related 
to the local wavevector by vlp = hkLP /ttilp^^ . 

For small wave vectors, \kx — kp\(^ <C 1, the Bogoliubov 
dispersion tends to a (lossy) sonic dispersion 

^Bog{kx)-uJp ^ Cs \kx-kp\^VLp{kx-kp)^i^- (11) 

with a sound speed Cs of fluid excitations such that 
"n^Lpc^ = J^gLP^LP^ and a Doppler shift term due to 
the background fluid flow at vlp- For large wavevectors, 
it recovers a (lossy) parabolic single-particle dispersion 

ZTUlp 

QLpriLP ^i^^^ (12) 

with a Hartree energy shift hriLpgLP due to the inter- 
actions with the condensate. Plots of the Bogoliubov 
dispersion at different spatial positions are shown in the 
small panels of Figs. |3] and [4j In the following of the 
paper, the will be used to interpret the numerical results 
for the evolution of excitation wavepackets and for the 
density-density correlations. 

Before proceeding, it is important to emphasize the 
conditions underlying the dispersion ([lO| of the collective 
excitations. In the coherent pump case (i), this form 
only holds at the end point of the upper branch of the 
hysteresis loop, while gaps and/or band sticking regions 



appear in all other cases because of the coherent pump 
locking the condensate phase^^ 

On the other hand, the presence of a soft excitation 
branch of sonic nature is guaranteed in the ballistic case 
(ii) by the fact that the condensate phase is completely 
free to evolve both in space and tim^^. However, as 
the density is spatially varying, the very concept of dis- 
persion relation holds only locally and is limited to ex- 
citations whose wavelength is small as compared to the 
characteristic length scale of the density profile. 



D. Truncated Wigner method 

A standard technique to go beyond the mean-field ap- 
proximation and include the fluctuations of the polari- 
ton field around its mean value is based on the Wigner 
representation of quantum fields. Within the truncated 
Wigner approximation, the dynamics of the quantum 
field problem can be described by a stochastic partial dif- 
ferential equation that can be numerically solved. This 
technique was first introduced in the context of quantum 
fluids in Ref.[44land its application to atomic condensates 
was fully developed and characterized in Ref. 27, Its ex- 
tension to the polariton case was developed in Ref. ^28l 
remarkably, the presence of loss and pump terms in the 
stochastic equations suppresses some of the difficulties of 
the truncated Wigner method and guarantees its accu- 
racy as long as interactions between individual polaritons 
are weak enough. If a spatial grid of real-space spacing 
Ax is used to numerically solve the stochastic partial dif- 
ferential equation, the weak interaction condition can be 
formulated as \g\ <^ j Ax. While this condition rules out 
the possibility of using the truncated Wigner method to 
study strongly correlated polariton states, e.g. in the po- 
lariton blockade regime^^ it does not hinder the study 
of quantum hydrodynamics effects such as the analog 
Hawking radiation, which originate from the collective 
dynamics of a large number of polaritons. 

The stochastic partial differential equations of the 
truncated Wigner method for the coupled exciton and 
polariton fields has the explicit form^^ 



ih 





hE{x,t) 



%(l^xp-i,) 



dt- 



V4 Ax 



where Ax is the spacing of the real-space grid and dW^ 
and dWc are complex valued, zero-mean, independent 
Gaussian noise terms, with white noise correlation in 
both space and time 

dW^{x,t) dWj{x',t) = 2Sa:,x' Sij dt (14) 

with z,j = {x, c}. In the numerical simulations, we re- 



construct the equilibrium Wigner distribution by first 
letting the system evolve to its steady state under the 
monochromatic pump, and then taking a large number 
of independent configurations by sampling the stochas- 
tic evolution at different times spaced by Tgam- In order 
to ensure statistical independence of the different real- 
izations, a large enough Tgam has to be taken such that 
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^sam7x,c ^ 1. In practice, a number of realization on the 
order of 10^ is used, taken at time intervals on the order 
of Tsam = IOTjTc- usual in Wigner approache d^^ * ^^ *, 
the stochastic averages over the configurations of differ- 
ent functions of the fields provide the expectation value 
of the corresponding symmetrically ordered operator. 

Inspired from previous numerical studies of acoustic 
black holes in atomic condensates^^, we shall look for the 
signature of analog Hawking radiation in the normalized, 
zero-delay correlation of the cavity field intensity defined 



.(2) 



(^00 ^ 00 ^ 



(15) 



Within the Wigner formalism, the different contributions 
to (15) have the form^^ 



1 



2 Ax 



(16) 



and 



(l + ^.,.O<l0c(x)P + |(/>c(xOP)^, (17) 



2 Ax 



where the {. . .)w averages indicate the classical averages 
over the different stochastic configurations of the field 
sampled at interval times Tgam- 



III. CREATING THE ACOUSTIC BLACK HOLE 

After having reviewed the technical tools to study the 
dynamics of the polariton quantum fiuid, we can now pro- 
ceed to discuss realistic configurations that can be used 
to generate an analog acoustic black hole with a large 
surface gravity. All simulations are performed by solv- 
ing the generalized Gross-Pitaevskii equation ^ using 
experimental parameters taken from the recent work in 
Ref.|33, 

The most straightforward configuration to generate an 
acoustic horizon was proposed by Solnyshkov et al.^^ and 
is illustrated in Fig.jsj^b): the presence of losses is respon- 
sible for a spatially decreasing density profile of a ballisti- 
cally fiowing polariton condensate and a correspondingly 
increasing velocity profile. Eventually, the polariton den- 
sity tends to zero and the fiow becomes necessarily super- 
sonic. The polariton density and speed in the pumped 
region is fixed by the pump beam parameters: if these 
are suitably chosen to have an initially sub-sonic fiow, a 
horizon must necessarily appear at some point. 

In spite of the simplicity of this configuration, some 
care has to be paid to the value of Hawking temperature 
that can be expected. For smooth horizons within the 



hydrodynamic limitP^, the Hawking temperature Th is 
determined by the surface gravity 



1 



d 



2cs{x) dx 



[vlp{x) - cl{x)] 



(18) 



according to Th = hn/kB^ with hp the Boltzmann con- 
stant. According to our numerical simulations and ana- 
lytical approximations'^, the characteristic length scale 
of the variation of the fiow parameters (density and 
flow speed) is inversely proportional to the loss rate 
^ (X VLp/^Lp. Direct combination of these two general 
facts shows that the achievable values of the Hawking 
temperature is limited from above by the loss rate: for 
the case in the figure, one indeed has fi/^ :^ 0.04 meV. 

A possible solution to overcome this difliculty, is to 
insert a narrow repulsive potential in the ballistic flow 
region, close to the edge of the pump spot. The efli- 
ciency of this technique is illustrated by the numerical 
calculations shown in Fig. [4]^b). Flow across the defect 
mostly occurs via tunneling processes through the poten- 
tial barrier: this produces a sudden drop of the conden- 
sate density and a corresponding sudden increase of the 
flow speed. As a result, a horizon appears in the vicinity 
of the defect center, separating a sub-sonic upstream re- 
gion from a super-sonic downstream one, with a surface 
gravity as high as hn 0^1.2 meV, i.e. about a factor ~ 30 
larger than the purely ballistic flow case of Fig. [3] 

Another difficulty of the purely ballistic flow conflgu- 
ration of Fig. [sjb) was the signiflcant spatial variation of 
the polariton density on both sides of the horizon, which 
may strongly distort the geometrical structure Hawking 
signal. The defect conflguration of Fig. [4]^b) appears to 
be favorable also in this respect: as the flow in the up- 
stream region is sub-sonic, no Bogoliubov-Cerenkov emis- 
sion can take place from the defect in this direction^^ and 
the polariton density remains almost flat across the whole 
flat-top region. The slow spatial variation of the density 
in the downstream super-sonic region is hardly avoided 
unless more complex laser conflgurations are used, but 
does not seem to prevent identiflcation of the Hawking 
effect. 

On the other hand, the use of a coherent pump raises 
the crucial issue of avoiding all those branch- stickin g and 
gap effects typical of non-equilibrium condensateP^^^ 
that would spoil the analogy with gravitational systems. 
As we have reviewed in Sec. II C[ choosing the value of 
the laser amplitude in the flat-top region in the vicinity 
of the end-point of the upper branch of the hysteresis 
loop is enough to have a linear, sonic dispersion of the 
low-/c dispersion of the form ([To| at all spatial positions 
around the horizoi]^. 

As it was illustrated in Fig. [2j keeping the system sta- 
bly on the upper branch of the hysteresis loop requires in- 
cluding a region of higher pump intensity upstream of the 
flat- top. Exception made for its role in switching on the 
optical bistability, the high pump intensity region play 
no role in the Hawking physics. As we have discussed in 
Sec. II C[ the Bogoliubov dispersion in the downstream 
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Figure 3: (Color online) (a) Scheme of the purely ballistic flow 
configuration to create a flowing polariton superfluid with an 
acoustic horizon, (b) Spatial profiles of the flow t'LP and the 
sound Cs speeds at steady state. The parameters and the spa- 
tial shape of the pump laser are the same as in Fig. [2] (c) Dis- 
persion of the collective excitations at two spatial positions: 
point Xp (left panel) is chosen to be in the sub-sonic flat-top 
region; point Xprobe (right panel) is located downstream of 
the horizon in the super-sonic region. Red solid (blue dashed) 
lines refer to positive (negative) norm Bogoliubov modes, (d) 
Spatial profile of the exciton density at a time t = 50 ps af- 
ter the arrival of the probe pulse. The probe laser parameters 
are chosen so to resonantly excite the Bogoliubov branch indi- 
cated as di^in in panel (c): frequency h{(jOs — cjp) = 0.13 meV, 
wavevector — kLp{xprohe) — — 1.55 /xm""*^. The probe spot 
has a gaussian shape centered at Xprobe = 450 /xm with a waist 
Ws — 20 /xm and a peak amplitude hF^ — 0.35 meV-/im~^/^. 
Its temporal shape is also gaussian of duration r = 10 ps. 
All numerical calculations are performed in a 800-/im simu- 
lation box with 2048 lattice points and absorbing boundary 
conditions at the edges of the box. 



Figure 4: (Color online) (a) Scheme of the ballistic flow 
configuration with a defect so to create an acoustic hori- 
zon with a large surface gravity. Panels (b-d) are analo- 
gous to the corresponding ones of Fig. [3] The parameters 
and the spatial shape of the pump laser are again the same 
as in Fig. [2] The defect consists of a square potential bar- 
rier of thickness Ldef = 2.5 /xm and is located at a position 
Xdei = 354 /xm close to the flat-top edge. The height of the 
defect potential is Vdef = 1 meV and acts on both exciton 
and photon components of the polariton field. The horizon 
appears in the close vicinity of the defect with a large sur- 
face gravity. The probe laser parameters are chosen so to 
resonantly excite the Bogoliubov branch indicated as (il,^n 
in panel (c): frequency h{uJs — ujp) — 0.5 meV, wavevector 
ks — khP (xprohe) = — 1.9 /xm""*^. The shape and the amplitude 
of the pump spot are the same as in Fig. [sjd). 



IV. WAVEPACKET SCATTERING 



region is guaranteed to be sonic at low wavevectors by 
the ballistic nature of the flow in this region. The flow 
speed is guaranteed to be super-sonic by the suppressed 
value of the density beyond the defect. 



The numerical simulations that we have presented in 
the previous section demonstrate the possibility to cre- 
ate configurations of a polariton flow through a horizon 
separating regions of sub- and super-sonic regimes with 
a sizable surface gravity. This configuration will be a 
powerful workbench for studies of Hawking phenomena. 
Inspired by related theoretical work on atomic conden- 
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sateJ^ and experiments in classical hydro dynamic^^, in 
this Section we shall first consider the scattering of a 
phonon wavepacket off the horizon, a process where the 
classical counterpart of the Hawking effect manifests it- 
self as an additional reflected wavepacket. This discus- 
sion will be the starting point for the next section, where 
we shall investigate the very Hawking effect, namely the 
conversion of zero-point quantum fluctuations into ob- 
servable radiation by the horizon. 

The scattering dynamics of Bogoliubov phonons hit- 
ting the horizon can be studied in terms of the generalized 
GPE far from the horizon, the propagation of weak 
wavepacket perturbations follows the Bogoliubov disper- 
sion discussed in Sec. |H C| albeit with spatial dependent 
flow and sound speeds. Figs. |3] to [6] illustrate the different 
aspects of this physics for the two previously mentioned 
cases, namely a smooth flow with a low surface gravity 
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Figure 5: (Color online) Snapshots of the spatial profile of the 
density modulation due to the propagating phonon wavepack- 
ets at different times (a) t = 50 ps, (b) t = 100 ps, and (c) 
t = 175 ps after the arrival of the probe pulse. Each panel 
displays the excitonic density modulation with respect to the 
time-independent steady state. The pump and probe param- 
eters are the same as for the purely ballistic flow shown in 
Fig.jsjd). The labels on the wavepackets refer to the Bogoli- 
ubov dispersions of Fig. [sj^c), while the arrows indicate the 
propagation direction of each wave packet. 
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Figure 6: (Color online) Snapshots of the spatial profile of the 
density modulation due to the propagating phonon wavepack- 
ets at different times (a) t = 50 ps, (b) t = 100 ps, and (c) 
t = 175 ps after the arrival of the probe pulse. Each panel 
displays the excitonic density modulation with respect to the 
time-independent steady state. The pump and probe parame- 
ters correspond to the defect configuration shown in Fig.[4|d). 
The labels on the wavepackets refer to the Bogoliubov disper- 
sions of Fig. [4]^c), while the arrows indicate the propagation 
direction of each wave packet. 



a defect configuration with a sizable surface gravity. For 
each configuration, the structure of the flow is illustrated 
in Figs.jsjb) and[4]^b), while the smaller (c,d) panels show 
the Bogoliubov dispersion at different spatial positions. 

Once the system is in its steady state, a low-frequency 
wavepacket perturbation is generated in the downstream 
super-sonic region by means of an extra coherent laser 
pulse with a spatio-temporal shape of the form 



-{X-Xs) /V 



i{ksX—ujst) 



(19) 



The probe wave vector (ks) and frequency {us) are cho- 
sen to be resonant with the Bogoliubov branch labeled 
as di^in in the (c) panels. The probe beam is centered at 
^probe; its waist Wg determines the spatial size of the gen- 
erated wavepacket (and its inverse width in /c-space). The 
temporal duration r has to be long enough not to excite 
the upper polariton branch nor the negative-norm Bogoli- 
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ubov branch. Its amplitude Fs is chosen weak enough to 
remain within the vahdity domain of a hnearized descrip- 
tion of excitations. 

The generated wavepacket shown in Figs, [sj^d) and 
[Ij^d) then propagates in the leftward direction against 
the horizon. A series of snapshots of this evolution are 
presented in Figs. |5] and [6j The difference between the 
behaviors of the weak and of a large surface gravity cases 
is apparent: in the former case, only two wavepacket s 
visibly emerge from the horizon at long times, located on 
the 

Uout ^iid di out Bogoliubov branches. On the other 
hand, in the latter case, an extra wavepacket appears on 
the d2^out branch. Remarkably, the Bogoliubov norm of 
this wavepacket is negative, which signals the occurrence 
of Hawking conversion processes. The identification of 
the different wavepackets in terms of their Bogoliubov 
branch is confirmed by the numerical measurement of 
their group velocity and of their wavevector. 

Of course, the physical reason why we did not observe 
a d2^out wavepacket in Fig. |5]is that the low value of the 
surface gravity strongly suppresses the amplitude of the 
Hawking process. As it was expected on general ground^ 
and then explicitly verifie d * the Hawking scattering 
amplitude for a Bogoliubov wavepacket of given carrier 
frequency uj scales down exponentially with the surface 
gravity according to the Boltzmann factor ex.-p{—27ruj/hz). 
On the other hand, excitation of the 6^2, ont branch would 
be totally forbidden if the fiow was everywhere sub-sonic 
and no horizon was present: in this case, in fact, energy 
could not be conserved in the Hawking scattering process. 



V. SPATIAL CORRELATIONS OF DENSITY 
FLUCTUATIONS 

The process studied in the previous section provides a 
solid evidence of the classical counterpart of the Hawk- 
ing effect, namely the possibility of inter-converting pos- 
itive and negative norm waves at the horizon. Now we 
shall proceed with the study of the Hawking effect stricto 
sensu, i.e. the emission of radiation by the horizon via 
the inter-conversion of zero-point quantum fiuctuations 
into observable radiation. 

To this purpose, one has to fully include in the model 
the quantum fiuctuations of the fields around the mean- 
field value predicted by the Gross-Pitaevskii equation. In 
this work, we adopt the truncated Wigner approach de- 
scribed in Sec. |IID[ where expectation values of quantum 
operators are calculated as classical averages over a suit- 
ably chosen stochastic partial differential equation. In- 
spired from previous theoretical^ and numericaP^l work, 
we shall consider the spatial correlation function of the 
quantum fiuctuations of the fiuid density: the smoking 
gun of Hawking radiation consist of a negative long- 
range correlation signal between points located on op- 
posite sides of the horizon. 

Numerical plots of the normalized density-density cor- 

(2) 

relation function Qc (x, x') — 1 are shown in the two pan- 




X - Xhor (Mm) 



Figure 7: (Color online) Color scale plot of the normalized 
spatial correlation function of photon density fluctuations 
gi^\x, x'). This quantity directly reflects in the intensity cor- 
relations of a near-field image of the emitted light from the 
microcavity. The upper (a) panel refers to the case of a purely 
ballistic flow considered in Figs, [s] and [s] The lower panel (b) 
refers to the case with a defect in the ballistic flow considered 
in Figs. [4] and [g] The colored dashed lines in (b) indicate the 
different Hawking features as discussed in the text. 



els of Fig. [7| for the two cases of a weak and strong sur- 
face gravity respectively. While in the upper panel, 
Fig. [7|^a) (weak surface gravity, hti ^ 0.04 meV) there 
is no clear feature emerging over noise, the lower panel. 
Fig. [7|^b) (strong surface gravity, fin ~ 1.2 meV) shows 
several features that can be interpreted a s a di rect con- 
sequence of the analog Hawking radiatioiP^l^ the two 
Bogoliubov excitations forming the Hawking pair are si- 
multaneously emitted from the horizon. At later times, 
the temporal correlation between the emission time of 
the Hawking partners reflects into a long-distance corre- 
lation in the density fiuctuations due to each member of 
the pair. While in the atomic case this correlation ex- 
tends at long times to indefinitely large distances from 
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the horizon, the finite decay rate of Bogohubov excita- 
tions in polariton fluids predicted in ([lO| restricts the 
correlation signal to a finite region of size v^/^lp around 
the horizon, being the group velocity of the Bogoh- 
ubov branch under consideration. Away from this region 
around the horizon, the Hawking phonons have decayed 
and the corresponding density fluctuations have disap- 
peared. The dashed lines in Fig. [7|^b) indicate the loca- 
tion at which one would expect the maximal density cor- 
relations due to the quantum vacuum emission process 
into the Uout - d2,out (white), the Uout - di^out (purple), 
the di^out — d2,out (yellow) pair of modes. 

The first process is the traditional Hawking emission 
process where the two quanta are emitted into oppo- 
site directions from the horizon. The main properties 
of the density correlation pattern due to this process is 
well captured by the non-dispersive quantum field the- 
ory on curved space-time as proposed in Ref. 51, For 
X — Xhor > and — Xhor < 0, the correlation is peaked 
along the half-line 



Xhor 



Xhor 



Vd - Cd 



(20) 



Here, Vd,u and Cd^u are the flow speed and speed of 
sound in the upstream and downstream regions, so that 
Vd — Cd and Vu — Cu are the group velocity of the 6^2, ont 
and Uout modes in the low-/c, sonic region. While in 
the atomic BEC case of Ref. T3" their value did not de- 
pend on space in the asymptotic regions away from the 
horizon, the driven-dissipative nature of polaritons is re- 
sponsible for the spatial dependence of the flow vlp and 
sound Cs speeds that is visible in Fig. [4]^b), mostly in 
the downstream region. However, as the space depen- 
dence remains moderate within the regions of thickness 
\vu - Cu\/^LP and |q - Vd\/^LP on either side of the 
horizon where correlations are significant, no appreciable 
curvature of the Hawking tongues is visible in Fig. [7|^b). 

The orientation of the purple and yellow dashed lines in 
Fig. [rj^b) is analogously obtained from the group velocity 
of the corresponding modes, namely 



hor 



■X, 



hor 



Vd^Cd 



and 



X 



Vd - Cd 



Vd + Cd 



(21) 



(22) 



The excellent agreement of the geometrical location 
of the numerically observed features with the analytical 
predictions (20-22) confirms our interpretation. Remark- 
ably^^, the peak value of the Uout — d2^out Hawking cor- 
relation signal is on the order of g^^^ - 1 0.92 x 10-^ 
which is not too far from the analytical prediction of 
Ref. HH As the surface gravity is here comparable to 
the interaction energy hgLpTiLP, the quantitative dis- 
crepancy can be easily traced back to the breakdown of 



the hydrodynamical approximation underlying the ana- 
lytical model^^. Thanks to the continuous- wave nature 
of the proposed experimental setting, the quantitatively 
low value signal to be observed can be overcome by a 
sufficiently long integration time. 

While agreement with the atomic case is good also for 
the Uout - di^out feature, the di^out - d2,out one has a 
different sign. A possible explanation of this behavior 
can be traced back to the non-universality of the Uout 
di^out back-scattering effect that is responsible for the 
conversion at the horizon of the standard Uout — d2^out 
Hawking correlation into the di^out — d2,out one. 

Before concluding, it is important to note that 

throughout the whole discussion, we have implicitly as- 

(2) 

sumed that the gc {x^x') correlation function is assumed 
to be measured at the same time: a rough estimate of the 
required temporal resolution of the detector is given by 
5t ~ 5x/cu4^ where Cu^ is the speed of sound in the po- 
lariton gas (on the order of 2/im/ps from Figs. |4| and 
5x is the spatial width of the correlation signal. Using 
the value (^x ~ 8 jiia taken from Fig. [7| one can esti- 
mate the needed temporal resolution to be on the order 
of 5t ~ 4ps, which is within the state-of-the-art of optical 
technology. 



VI. CONCLUSIONS 

In this article we have presented a comprehensive study 
of classical and quantum hydrodynamic properties of 
acoustic black holes in superfluids of exciton-polaritons in 
semiconductor microcavities. Inspired from on- going ex- 
perimental research, we have identified and characterized 
one-dimensional polariton wire devices as model systems, 
where the polaritons injected by a single monochromatic 
laser beam with a suitable spot profile generate a flow 
configuration showing an acoustic black hole horizon. In- 
serting a repulsive potential defect provides a dramatic 
enhancement of the analog surface gravity at the hori- 
zon. Even if our calculations have been performed for 
the specific case of a semiconductor microcavity in the 
strong-coupling regime, all conclusions straightforwardly 
extend to a generic planar cavity filled by a nonlinear 
optical medium, as originally suggested in Ref. |20l or 
photonic crystal polariton states with engineered mode 
dispersion^^ 

In analogy with previous works on acoustic black holes 
in atomic condensates and classical hydrodynamics of 
surface waves, two experiments have been proposed and 
numerically simulated to assess the efficiency of Hawking 
wave-conversion processes at the horizon. When a coher- 
ent wavepacket of Bogoliubov phonons is incident on the 
horizon, the Hawking effect is visible as the production 
of an extra emerging wavepacket on the negative-norm 
branch. This is the classical evidence for the Hawk- 
ing process taking place at the horizon. On the other 
hand, the usual Hawking radiation is a purely quantum 
effect consisting of the conversion of zero-point fluctua- 
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tions into a stream of observable phonons at the horizon: 
the correlation between the emission times of a phonon 
and its Hawking partner is visible as a long-distance cor- 
relation between the density fluctuations on either side of 
the horizon. The density correlations provide the clear- 
est experimental signature of Hawking radiation in su- 
perfluids: in the present polariton case, in-cavity density 
correlations directly transfer into the secondary emission 
from the cavity photon, and can be observed as analogous 
correlations in the intensity fluctuations of the near-field 
emission. As soon as a sufficient signal over noise ratio 
is achieved, an experiment along the proposed lines will 
hopefully provide a clear evidence of a fundamental ef- 
fect of quantum field theory that so far has eluded all 
experimental observation, overcoming the barrier to di- 



rect observation imposed by the negligibly small Hawking 
temperature of astrophysical black-holes. 
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Restricting to the bottom of the lower polariton branch, 
it is easy to see that the polariton density can be ex- 
pressed as riLP — |0x|^ + 10c 1^, and the polariton speed 
is defined as vlp = jLp/riLP in terms of the polariton 
current Jlp = 4^1'^LP^'ipLP - ^Plp^^Ip] ^ jc = 

2im VV^c — ipc VV^*]. Given its very heavy mass, ^ 
mc, the exciton fraction does not significantly contribute 
to the polariton current. 

In practical simulations, we keep the pump amplitude at a 
value slightly above the turning point to allow a faster and 
more robust convergence towards the steady state. Even 
if the resulting Bogoliubov dispersion shows a small gap, 
its amplitude is however much smaller than the polariton 
broadening jlp, and hence has no influence on the final re- 
sults, which further supports the robustness of our scheme. 
While the qualitative structure of the Hawking correlation 
pattern appears quite similar, quantitative comparison of 
our results with the ones in Ref. [24] is hardly made, as 
the scale of the correlation signal is missing there. Differ- 
ently from that paper, the negative sign of the correlation 
signal reported in the present work matches the one ex- 
pected from theoretical and numerical work in equilibrium 
condensatePI^. This sign difference is most likely due to 
the fact that in Ref. '2^ fluctuations are induced by some 
classical disorder present in the super-sonic region. 



